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HD 要 : 土壤 电导 率 (Electrical conduetivity ,BC) 是 评价 土壤 盐 渍 化 的 重要 指标 。 通 过 实测 新 疆 艾 
比 湖 湿地 自然 保护 区 土壤 EC 及 可 见 光一 近 红 外 光谱 数据 ,利用 波谱 响应 技术 模拟 Landsat 8 OLI, 
Sentinel 2 Sentinel 3 卫星 的 宽 波段 数据 。 构 建 宽 波段 模拟 数据 及 其 5 种 预 处 理 后 的 三 维 光 谱 指 数 
( Three-dimensional spectral index, TDS/) ,采用 梯度 提升 回归 树 算法 (Gradient boosting regression 
tree, GBRT) 建立 3 种 卫星 土壤 EC 估算 模型 HLM MATOS 后 模型 精度 的 变化 。 结 果 表 明 : 在 
不 同 土壤 EC 条 件 下 ,3 种 卫星 具有 相似 的 光谱 趋势 , 均 在 红 、 近 红外 波段 附近 反射 率 较 高 ;TDSIJ 与 
土壤 EC 相关 性 基本 均 在 0.4 以 上 ,最 大 程度 保留 了 与 土壤 EC 敏感 度 高 的 红 、 绿 、 蓝 、 近 红外 、 短 波 
红外 波段 信息 ;GBRT 对 于 土壤 EC 估算 能 力 表 现 突出 ,3 种 卫星 对 土壤 EC 的 最 佳 预测 精度 R 分 
别 为 0.831、0.847、0.903, 在 加 入 TDSI 后 ,R 分 别提 高 至 0.835 .0.857 .0.935 ,综合 分 析 发 现 ,Sen- 
tinel 3 对 土壤 EC 估算 效果 最 佳 (R =0. 935 , 均 方 根 误差 RMSE =2.986 mS . em , 赤 池 信息 准则 
AIC =57. 500)。 通 过 利用 波谱 响应 技术 结合 TDS 深度 挖掘 波段 间 的 协同 信息 ,采用 GBRT 验证 了 


量 监测 与 防 控 提供 有 利 指导 。 
关键 词 : 
文章 编号 : 


土壤 盐 渍 化 是 中 亚 乃 至 全 球 干旱 区 重要 的 环境 
问题 之 一 , 它 不 仅 阻碍 植物 的 正常 生长 ,影响 农业 系 
统 的 生产 力 , 同时 还 在 一 定 程度 上 造成 土地 退 
4t! 7 。 新 疆 维吾尔 自治 区 位 于 我 国 西北 干旱 区 ， 
是 国家 重要 的 粮食 产业 基地 ,为 确保 稳定 的 粮食 产 
出 和 粮食 安全 ,及 时 监测 和 防 控 土 壤 盐 渍 化 显得 十 
分 必要 。 因 此 ,实现 干旱 区 土壤 盐分 的 精准 监测 
和 有 效 管理 对 该 区 域 可 持续 发 展 和 维持 生态 稳定 具 
有 重要 意义 。 

地 物 光 谱 特征 作为 遥感 搁 术 及 应 用 中 的 重要 物 
理 基础 与 理论 依据 ,是 反映 地 物 信息 的 有 效 手段 ,已 
有 众多 国内 外 学 者 在 土壤 植被 .水体 等 领域 开展 地 
物 光谱 相关 研究 工作 ““ 。 地 物 光 谱 能 够 准确 反映 
地 物 的 光谱 特性 ,常见 的 光谱 研究 形式 主要 有 光谱 
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不 同 卫星 对 土壤 EC 的 估算 效果 ,二 者 相 结合 可 以 有 效 提 升 模型 预测 精度 ,为 干旱 区 土壤 盐 渍 化 定 


波谱 响应 ; 土壤 电导 率 ; 光谱 预 处理 ; 三 维 光 谱 指数 ; 梯度 提升 回归 树 
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预 处 理 法 ,光谱 指数 法 等 ,光谱 预 处 理 能 够 有 效 减 少 
背景 噪声 ,增强 光谱 特征 ,而 通过 构建 光谱 指数 则 可 
以 深度 挖掘 光谱 间 的 协同 信息 ` 减 少 信息 元 余 ,并 在 
一 定 程度 上 提升 对 地 物 的 响应 程度 '"。 但 现 阶段 
光谱 指数 大 多 仅 考 虑 两 波段 间 的 相互 关系 ,这 可 能 
会 忽略 部 分 地 物 敏感 信息 ,造成 估算 偏差 ,而 三 维 光 
谱 指 数 在 增加 一 维 波段 的 基础 上 ,增加 了 光谱 有 效 
言 息 , 尽 可 能 地 保证 了 估算 精度 ,是 光谱 指数 研究 的 
新 思路 。 

传统 的 土壤 盐分 监测 依靠 野外 实地 采样 实验 
室 室 内 理化 分 析 的 方法 ,该 方法 虽 具 有 和 较 高 测量 精 
BE ,但 难以 满足 大 太 度 上 的 监测 需求 。 遥 感 技术 依 
靠 其 覆盖 范围 广 .时 间 连 续 性 强 等 优点 ,成 功 实现 宏 
观 .多 时 相 的 土壤 盐分 监测 ,为 土壤 盐 渍 化 的 监测 和 
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防 控 提供 了 有 利 保障 '"。 当 前 ,相关 研究 主要 集中 
于 利用 不 同 遥感 数据 针对 不 同 区 域 的 土壤 盐 渍 化 信 
息 反 演 ,并 将 实测 样 点 数据 与 遥感 反 演 数据 进行 精 
度 评 佑 ,结合 不 同 模型 实现 土壤 盐 渍 化 信息 监 
3^. Rep ,部 分 学 者 还 探讨 不 同 算法 对 土壤 盐 
分 的 估算 精度 ,进而 选择 最 佳 算 法 实现 土壤 盐分 的 
精确 估算 777 。 在 以 上 研究 分 析 中 发 现 , 利 用 遥感 
数据 反 演 土壤 盐 省 化 信息 是 切实 可 行 的 监测 手段 ， 
且 研 究 结果 表明 ,土壤 电导 率 与 土壤 盐 溃 化 具有 较 
高 的 相关 性 ,可 在 一 定 程 度 上 作为 衡量 盐 渍 化 程度 
的 指标 "| 。 通 过 既得 文献 梳理 归纳 发 现 ,利用 地 物 
实测 光谱 数据 模拟 不 同 卫星 传感器 波段 信息 ,将近 
感 与 遥感 技术 协同 融合 ,进而 实现 土壤 盐分 定量 监 
测 的 研究 尚 不 多 见 。 新 疆 地 区 土壤 盐 溃 化 类 型 丰 
富 ,程度 复杂 ,在 全 球 干 旱 区 的 土壤 盐 渍 化 研究 中 非 
党 典型 ,加 之 该 区 域 下 垫 面 均 一 ,数据 获取 质量 较 
高 ,因此 ,尝试 利用 地 物 实测 光谱 与 遥感 卫星 传 感 融 
进行 波谱 响应 研究 ,以 期 为 全 球 干旱 区 的 土壤 盐 渍 
化 精确 监测 与 数字 制图 提供 依据 。 

本 研究 以 新 疆 维 吾 尔 自治 区 艾 比 湖 湿地 自然 保 
护 区 为 研究 靶 区 ,通过 实测 光谱 数据 模拟 Landsat 8 
OLI „Sentinel 2 „Sentinel 3 宽 波 段 ,利用 多 种 预 处 理 
方式 构建 三 维 光 谱 指数 (Three-dimensional spectral 
index, TDSI) , 六 选 最 优 光 谱 指 数 。 将 实测 土壤 电 
导 率 (Electrical conductivity , EC ) 数据 结合 梯度 提升 
回归 树 算 法 ( Gradient boosting regression tree, 
GBRT) 建 立 土壤 EC 估算 模型 ,并 比较 三 种 卫星 对 
土壤 EC 的 估算 潜力 ,为 干旱 区 土壤 盐 渍 化 监测 与 
防 控 提供 有 利 指导 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

艾 比 湖 湿地 地 处 博 尔 塔 拉 蒙古 自治 州 ,位 于 天 
山北 蔓 准 葛 尔 盆地 西南 隅 43°38" ~ 45°52'N, 
79°53" ~ 85°02'E) ,是 典型 的 干旱 区 尾 闻 盐湖 湿地 。 
该 区 域 全 年 平均 降水 量 约 为 105 mm ,年 均 潜在 蒸发 
量 1 315 mm, 年 均 温 8.3 °C ,土壤 类 型 以 荒漠 土 IK 
漠 士 为 主要 , 主要 植被 有 树 柳 (Tamarix) ik AK 
( Halocnemum strabliaceum) ,$h Jf. ( Kalidium folia- 
um) 等 “”。 由 于 处 在 大 陆 腹地 ,加 之 特殊 的 地 
EWM EEKMA WE ERR ABR. SA 


然 和 人 为 因素 的 影响 , 艾 比 湖 湖面 面积 持续 缩减 ,十 
壤 盐 渍 化 程度 不 断 加 剧 ,严重 危害 当地 的 生态 环境 
与 农业 生产 。 
12 土壤 样品 采集 与 土壤 EC 分 析 

2018 年 7 月 研究 团队 依据 研究 区 景观 类 型 , 选 
取 绿 洲 、 充 漠 以 及 交错 带 等 景观 ,按照 五 点 采样 法 采 
集 0 ~10 cm 土壤 样品 共计 40 个 ,用 样品 袋 封装 带 
回 实验 室 。 经 过 自然 干燥 后 ,研磨 通过 2 mm FL, 
并 将 每 个 样品 分 为 2 份 ,分 别 用 来 测定 土壤 EC 和 
可 见 光 一 近 红 外 光谱 数据 (Visible and near infrared , 
VIS-NIR) 。 过 得 后 的 土壤 样品 按照 土 水 质量 比 1:5 
的 比例 提取 土壤 浸出 液 ,使 用 德国 Wissenschaftlich 
Technische Werkstitten 公司 生产 的 Cond 7310 型 十 
壤 测 试 仪 测定 土壤 EC。 
1.3 光谱 采集 

利用 美国 Analytical Spectral Devices 公司 生产 
的 ASD FieldSpec 3 型 光谱 仪 (波段 350 ~2 500 nm) 
在 暗室 内 测定 样品 的 VIS-NIR 数据 。 将 制备 好 的 
40 个 土壤 样品 分 别 装 和 直径 为 12 cm ZR 3 cm 盛 样 
MP ,将 表面 刮 平 。 光 源 为 50 W 的 讽 素 灯 , 距 土壤 
样品 表面 50 cm ,天 顶 角 为 153“。, 探 头 至 待 测 样品 的 
表面 距离 为 10 cm, 每 次 光谱 测定 之 前 均 进 行 白 板 
标定 以 消除 暗 电 流 的 影响 "。 本 实验 各 土壤 样品 
于 8 个 方向 共 采 集 24 条 光谱 曲线 , 取 24 条 光谱 曲 
线 的 算术 平均 值 作为 该 土 样 的 原始 光谱 数据 。 
14 卫星 宽 波 段 光 谱 模 拟 及 预 处 理 

使 用 ENVI5.5 软件 中 内 置 的 光谱 重 采 样 功能 ， 
利用 Landsat 8 OLI, Sentinel 2 „Sentinel 3 卫星 各 自 的 
波谱 响应 函数 (图 1) ,将 实测 VIS-NIR 数据 模拟 重 
采样 为 每 个 卫星 传感器 的 宽 波段 反射 率 数 据 , 再 利 
用 ENVI 5. 5, Matlab R2016b 软件 对 其 进行 一 阶 微 
分 (First derivative, FD) 二 阶 微分 (Second deriva- 
tive, SD) .连续 统 去 除 (Continuum removal , CR ) 、 吸 
光度 变换 (Absorbance, ABS) 多 元 散射 校正 ( Mult- 
ivariate scattering correction, MSC)5 种 预 处 理 。 光 谱 
模拟 计算 公式 为 : 


(1) 


式 中 :有 为 模拟 卫星 宽 波段 的 反射 率 ;Aw Am 为 传 
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图 1 3 种 卫星 波谱 响应 函数 图 


Fig.1 Spectrum response function of three satellites 


感 器 光谱 探测 的 起 始 和 终止 波长 ;S, 为 传感器 在 和 
波长 处 的 波谱 响应 函数 值 ;R, 为 土壤 光谱 在 和 波长 
处 的 反射 率 。 
1.5 三 维 相关 系数 

光谱 指数 通过 波段 间 简 单 组 合 即 可 提高 光谱 信 
息 的 利用 率 ,降低 信息 宛 余 度 ,最 大 程度 上 保留 土壤 
EC 的 敏感 光谱 信息 ""“] 。 光 谱 指数 与 土壤 EC 间 
的 三 维 相关 系数 可 以 提供 关于 土壤 EC 的 3 个 独立 
波段 组 合 的 综合 信息 ,进一步 考虑 波段 间 的 协同 作 
用 ,增加 相关 信息 量 ,为 后 续 提 升 模型 预测 精度 奠定 
基础 。 本 研究 构建 两 种 TDSI ; 


TDSI, = (R, Rah (Ry -R,;) (2) 


TDSI, = (Ry RA [Ru 一 Ri ) = 
(Ry - R4) ] (3) 


SUP Ry, Ry Ra 分 别 为 波长 在 AAAs 处 的 反 
射 率 。 

光谱 指数 与 土壤 EC 的 相关 性 计算 以 及 三 维 相 
关系 数 图 的 呈现 均 在 Matlab R2016b 中 实现 。 
1.6 集合 划分 与 建 模 方法 

本 研究 采用 基于 联合 x-y 距离 方法 (Simple set 
portioning based on joint x-y distance ,SPXY ) 对 40 个 
样本 进行 划分 ,该 方法 能 够 有 效 履 盖 到 所 有 样本 , 增 
加 模型 预测 稳定 性 。 基 于 SPXY 方法 是 综合 考量 了 
光谱 反射 率 向 量 与 土壤 EC 向 量 的 情况 下 ,确保 样 
本 在 上 述 两 个 向 量 空间 中 具有 相同 的 权重 ,本 文 最 
终 确定 了 30 个 建 模 集 样本 ,10 个 验证 集 样本 。 


型 ,该 算法 是 由 FRIEDMAN 提出 的 一 种 迭代 的 决策 
MRO ,由 在 通过 集成 若干 个 弱 学 习 器 ,经 过 多 
次 迭代 最 终 组 合 而 成 一 个 强 学 习 器 ,每 一 次 迭代 都 
在 上 一 次 迭代 的 基础 上 减少 残 差 ,并 改进 上 一 次 的 
计算 结果 ,同时 根据 残 差 减少 的 梯度 方向 来 构建 新 
的 组 合 模型 。 模 型 预测 的 运算 机 制 大 体 为 :输入 一 
个 样本 实例 ,首先 对 其 赋予 一 个 初始 值 ,之 后 遍历 每 
一 棵 决策 树 ,每 棵 树 都 会 有 不 同 的 权 值 , 并 且 对 预测 
值 进行 相应 地 调整 修正 ,最 终 的 结果 是 将 每 一 棵 决 
策 树 的 结果 累加 得 到 最 后 的 预测 结果 ,计算 公式 如 
F: 


f(x) = Fy +B,T,(X) + B,T,(X) 9 + 
BuTu(X) (4) 


SUH ey AIME IB 为 权重 系数 ;7T(X) 为 每 一 棵 决策 
树 的 结果 。 
1.7 建 模 精 度 评价 

为 量化 模型 的 拟 合 度 ,评价 模型 的 估算 能 力 ,本 
人 研究 选用 决定 系数 (Coefficient of determination, 
R^) . 均 方 根 误差 ( Root mean squared error , RMSE ) 和 
赤 池 信息 准则 (Akaike Information Criterion , AIC ) 3 
种 参数 来 进行 评估 。47C 准则 是 检验 统计 模型 拟 合 
度 的 一 项 指标 ,其 基于 入 的 概念 ,平衡 模型 构建 中 的 
偏差 和 方差 ,一般 认 为 ,47C 值 越 低 代表 模型 拟 合 效 
spp". AIC 计算 公式 如 下 : 


SSR 


AIC = nlog (> ) + 2k (5) 


AP in 是 样本 量 ;SSR 是 模型 拟 合 残 差 平 方 和 ;k 是 


本 研究 采用 GBRT 算法 构建 土壤 EC 估算 模 


模型 自由 变量 个 数 。 
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2 BRAM 


2.1 样本 土壤 EC 统计 分 析 

表 1 为 土壤 EC 的 特征 统计 描述 。 全 样本 土壤 
EC 均值 为 7.240 mS + cm” ,标准 差 为 7. 586 , 建 模 
集 与 验证 集 各 自 均 值 分 别 为 7.175 mS + em, 
7.860 mS * cm- ,各自 标准 差分 别 为 7.270、8.791， 
全 样本 均值 与 标准 差 均 介 于 建 模 集 与 验证 集 之 间 。 
2.2 不 同 卫 星 宽 波段 模拟 反射 率 分 析 

图 2 为 不 同 土壤 EC 对 3 种 卫星 宽 波 段 模拟 反 
射 率 的 影响 。 分 析 图 2 可 以 得 知 , 随 着 土壤 EC 的 
增加 ,Landsat 8 OLI、Sentinel 2、Sentinel 3 各 波段 的 
反射 率 也 相应 增高 ,其 中 反射 率 最 低 波段 均 为 BI 
(0.433 ~ 0. 453 hm 、0. 429 ~ 0. 457 hm、0. 391 ~ 
0.409 jum) ,而 反射 率 最 高 波段 分 别 为 B6 (1.561 ~ 
1.651 jua), Bil (1. 562 ~ 1. 663 am), B21 
(0.997 ~1.017 um) 3 种 卫星 光谱 反射 率 整体 上 
均 呈 升 高 趋势 ,在 Landsat 8 OLI 与 Sentinel 2 中 , i 
波 红外 波段 (B7、B12 ) 反 射 率 均 略微 下 降 ,Sentinel 
3 反射 率 则 持续 升 高 。3 种 卫星 均 存 在 反射 率 陡 增 
的 波段 ,Landsat 8 OLI 为 B4 (0. 630 ~0.680 um), 
B5 (0.845 ~0.885 jun) ,Sentinel 2 为 B4 (0.641 ~ 
0.690 um) ,B10 (1.355 ~1. 393 um) ,Sentinel 3 为 
B7 (0.613 ~0.628 um) ,B12 (0.748 ~0.760 um), 


表 1 土壤 样品 EC 特征 统计 /mS .cm 
Tab.1 Statistical characteristics of soil electrical 


conductivity of soil samples / mS - cm! 


B17 (0.853 ~0.873 hm)。 通 过 分 析 发 现 ,不 同 卫 
星 反 射 率 最 高 . 陡 增 波段 大 致 均 位 于 红 光 、 近 红外 附 
近 , 这 与 前 人 研究 结果 一 致 -”| ,主要 是 由 于 土壤 
中 OH” \Si-OH 键 以 及 土壤 硅 酸 盐 矿 物 及 其 他 阳 离 
子 的 影响 。 
2.3 三 维 相 关系 数 分 析 

在 Matlab R2016b 软件 中 ,计算 3 种 卫星 基于 
宽 波段 模拟 数据 及 5 种 预 处 理 方式 与 土壤 EC 的 相 
关系 数 , 最 终 得 到 6 种 光谱 形式 下 的 TDSI, .TDSI, 
的 三 维 相 关系 数 。 

通过 表 2 可 知 , 在 3 种 卫星 不 同 光谱 形式 的 三 
维 相 关系 数 中 , 仅 基于 CR 预 处 理 下 Landsat 8 OLI 
卫星 的 TDSI, , TDSL, 以 及 Sentinel 2 卫星 的 TDSI, 
在 0.01 水 平 上 不 显著 (显著 性 检验 立 值 P,, = x 
0.403)。 在 Landsat 8 OLI 卫星 中 ,基于 ABS 预 处 理 
下 三 维 相 关系 数 最 高 ,分 别 为 0. 625, - 0. 626; 在 
Sentinel 2 Sentinel 3 卫星 中 ,基于 SD 预 处 理 下 TDSI, 
三 维 相关 系数 均 最 高 ,分 别 为 -0.634、- 0.720 , 基 
于 ABS 预 处 理 下 7DS 三 维 相 关系 数 同样 最 高 ,分 
别 为 -0.600、-0.678。 已 有 研究 表明 ,土壤 盐分 在 
蓝 红 光 、 近 红外 、 短 波 红外 附近 具有 较 高 的 敏感 程 
pe) ,而 通过 分 析 表 2 中 通过 指数 所 条 选 出 的 波段 
可 以 发 现 , 大 部 分 波段 同样 分 布 在 上 述 范 围 内 ,少许 
分 布 在 绿 光 附近 ,这 与 已 有 研究 结果 相 吻 合 ,充分 映 
证 了 三 维 光 谱 指 数 在 兼顾 遥感 光谱 机 理 的 同时 ,最 
大 程度 上 保留 了 土壤 EC 的 相关 信息 ,具有 一 定 的 
合理 性 与 科学 性 。 图 3 列 出 3 种 卫星 代表 性 的 三 维 
相关 系数 图 。 


样本 类 型 ”样本 数 ”最 大 值 ”最 小 值 均值 标准 差 " 
全 集 40 28.079 0.059 7.240 7.586 2.4 不 同 卫星 土壤 EC 建 模 分 析 
建 模 集 30 28.079 0.059 7.175 7.270 将 3 种 卫星 宽 波 段 模拟 数据 及 预 处 理 数据 分 别 
验证 集 10 27.309 0.466 7.860 ^ 8.791 作为 土壤 EC 估算 模型 的 独立 变量 ,实测 土壤 EC 数 
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图 2 
Fig.2 


波段 号 
不 同 土壤 EC 对 卫星 模拟 波段 的 响应 


Response of different soil EC to satellite simulated bands 


波段 号 
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Tab.2 Three-dimensional correlation parameters of three satellites under different pretreatment 


卫星 光谱 形式 TDSI, 相关 系数 (7) TDSI, 相关 系数 (7) 
Landsat 8 OLI R B1,B3,B2 -0.594 B2 B3 B1 -0.595 
FD B1 ,B6 B2 一 0. 599 B1 、B2 .B6 一 0. 604 
SD B2 \B7 \BS 0.423 B6 \B3 、B4 -0.415 
CR Bl B1 B1 0 BI B1 ,Bl 0 
ABS B2,B3,Bl 0.625 B2,B3,Bl — 0.626 
MSC B1 B3 ,B2 一 0. 594 BI,B3,B2 0.594 
Sentinel 2 R B1 B3 B2 -0.567 B1 B3 B2 0.567 
FD B1 .B10 ,B2 -0.595 B1 .B10 ,B2 0.589 
SD B2 , B6 ,B4 -0.634 BA ,B8A ,B3 -0.583 
CR B3 .B4 B2 - 0.363 B5 .B8 .B2 0.518 
ABS B2 .B1 ,B3 一 0. 599 B2,B3,Bl — 0.600 
MSC B1 ,B2,B3 —0.567 BlI,B3,B2 0.567 
Sentinel 3 R B3 .B5 .B4 -0.675 B3 B5 .B4 0.673 
FD B17 ,B18 ,B8 -0. 583 B8 ,B18,B9 -0.597 
SD B5 ,B10,B7 -0.720 B10 \B8 ,B12 0.579 
CR B4 .BS \B2 一 0. 468 B5 、B17 ,B12 -0.565 
ABS B3 .BS „B4 -0.684 B4 .BS ,B3 -0.678 
MSC B3 B5 .B4 -0.675 B3 B5 „B4 0.673 
0.720 
0.625 
0.432 
1 0.375 
: 0.144 
n 3 0.125 : 
ax 7 
i 2 -0.125 -0.144 
7 
Te -0.375 -0.432 
5 6 
Be 43 3 4 5 
$ 2112 ^ 5 -0.625 -0.720 
TDSI, Sentinel 2 SD 
0.60 0.678 
. 0.36 0.407 
: 0.12 0.136 
0125 oP 6 : i 
R 8 
-0.125 #10 -0.12 -0.136 
12 
12 -0.36 -0.407 
10 12 
dk 86 810 0.678 
4 4 6 -0.60 -0. 
2 
Be 2 $97 


TDSI, Landsat 8 OLI ABS 


TDSI, Sentinel 2 ABS 
图 3 土壤 EC 与 TDSI 的 三 维 相关 系数 图 


Fig.3 Three-dimensional correlation coefficient of soil EC and three-dimensional spectral index 


TDSI, Sentinel 3_ABS 


据 作 为 响应 变量 ,结合 GBRT 算法 构建 土壤 EC fi 
算 模型 ,分析 3 种 卫星 不 同 光 谱 形式 对 土壤 EC 的 
估算 效果 ,模型 结果 、 精 度 评价 如 表 3 所 示 。 


由 表 3 可 知 ,Landsat 8 OLI 卫星 基于 FD 预 处 理 
的 土壤 EC 估算 精度 最 佳 (R^ = 0. 831, RMSE = 
4.484 mS * em ^! , AIC =29.034) ,基于 CR 预 处 理 的 
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表 3 3 种 卫星 宽 波 段 模拟 土壤 EC 估算 结果 


Tab.3 Results of soil EC estimation for broadband simulation of three satellites 


建 模 集 验证 集 
卫星 光谱 形式 参量 
RMSE / mS * em ^! R RMSE / mS * em ^! R? AIC 
Landsat 8 OLI R 7 4.422 0.774 6.627 0. 434 32.42] 
FD 7 3.884 0.829 4.484 0.831 29.034 
SD 7 3.570 0. 823 5.198 0.772 30.317 
CR 7 5.478 0.521 6.939 0. 320 32.825 
ABS 7 4.424 0.781 6.661 0. 431 32.471 
MSC 7 3.796 0. 794 4.775 0.792 29.579 
Sentinel 2 R 13 4.204 0.822 6.605 0. 436 44.398 
FD 13 3.662 0.855 5.084 0. 824 42.124 
SD 13 3.356 0. 864 4.931 0. 847 40.989 
CR 13 4.197 0.784 5.780 0. 697 43.239 
ABS 13 4.202 0. 82 6.493 0. 461 44.248 
MSC 13 3.556 0.876 4.461 0.810 41.858 
Sentinel 3 R 21 4.771 0.703 6.900 0.359 60.777 
FD 21 3.501 0. 822 4.650 0. 823 57.349 
SD 21 2.263 0.946 3.515 0. 903 54.918 
CR 21 3. 405 0. 866 5.239 0. 686 58.385 
ABS 21 4. 891 0. 693 6.779 0. 399 60. 623 
MSC 21 3.516 0. 837 4.258 0. 832 56.584 


IER 为 决定 系数 ,RMSE 为 均 方 根 误差 ,A1C 为 赤 池 信息 准则 。 下 同 


土壤 EC 估算 精度 最 低 ( 尼 =0. 320, RMSE =6. 939 
mS .cm ,AIC =32. 825) ,各 光谱 形式 的 估算 精度 
依次 为 FD > MSC > SD >R > ABS > CR, Sentinel 2 
卫星 基于 SD 预 处 理 的 土壤 EC 估算 精度 最 佳 (R = 
0. 847, RMSE 24.931 mS .cm , AIC = 40. 989) , 
基于 R 的 土壤 EC 估算 精度 最 低 (R = 0. 436, 
RMSE =6.605 mS * cm ^! ,AIC = 44. 398) ,各 光谱 形 
式 的 估算 精度 依次 为 SD > FD > MSC > CR > ABS > 
R, Sentinel 3 卫星 基于 SD 预 处 理 的 土壤 EC 估算 
精度 最 佳 (R* =0. 903, RMSE 23.515 mS * em ^! , 
AIC 254.918) ,基于 R 的 土壤 EC 估算 精度 最 低 
(R? = 0.359, RMSE = 6. 900 mS + cem ^! , AIC = 
60.777) ,各 光谱 形式 的 估算 精度 依次 为 SD > MSC 
>FD > CR > ABS > R。 综 合 比 较 最 佳 估 算 精 度 后 认 
为 ,Sentinel 3 对 于 土壤 EC 的 估算 效果 最 好 ,其 次 为 
Sentinel 2 ,最 后 为 Landsat 8 OLI。 图 4 为 3 种 卫星 
在 构建 土壤 EC 估算 模型 中 各 波段 的 重要 性 ,分 析 
发 现 3 种 卫星 重要 性 较 大 的 波段 主要 是 蓝 波 段 、 红 
波段 . 近 红 外 波段 .短波 红外 波段 以 及 绿 波段 ,这 与 
前 文 所 得 发 现 一 致 。 

TDSI 能 够 提供 与 土壤 EC 更 敏感 的 相关 信息 ， 


本 研究 为 减少 建 模 过 程 中 敏感 度 较 低 波段 对 建 模 精 
度 的 影响 ,将 各 光谱 形式 下 通过 显著 性 检验 的 三 维 
光谱 指数 作为 光谱 参量 加 入 建 模 当中 ,分 析 模 型 效 
TR ,精度 评价 如 表 4 所 示 。 

在 加 入 TDSI 后 构建 的 3 种 卫星 不 同 光谱 形式 
的 土壤 EC 估算 模型 ,R 较 之 前 均 有 所 提高 , 且 均 达 
到 0.6 以 上 。 综合 分 析 发 现 ,3 种 卫星 基于 R 的 土 
HE EC 估算 模型 精度 提升 幅度 最 为 明显 ,其 中 Senti- 
nel 2 卫星 精度 提升 82. 3496 ,提升 幅度 最 大 ,其 次 为 
Sentinel 3 卫星 , #25+ 76.88% ,最 后 为 Landsat 8 OLI 
卫星 ,提升 72.81% ;基于 ABS 预 处 理 的 模型 提升 幅 
度 也 较 大 ,分 别提 升 66. 59% ,74. 44% 、68. 21% 。 
Landsat 8 OLI 卫星 基于 FD 预 处 理 的 土壤 EC 估算 
精度 最 佳 ( 尼 = 0. 835, RMSE 24.348 mS ，cm 1， 
AIC 232.765) ,基于 ABS 预 处 理 的 土壤 EC 估算 精 
度 最 低 (R = 0. 725, RMSE = 5. 601 mS * em^!, 
AIC =34.965) ,各 光谱 形式 的 估算 精度 依次 为 FD > 
MSC > SD >R > ABS, Sentinel 2 卫星 基于 SD 预 处 
理 的 土壤 EC 估算 精度 最 佳 ( 尼 = 0. 857, RMSE = 
4.596 mS + cem ^! , AIC =45.247) ,基于 CR 的 土壤 EC 
估算 精 度 最 低 ( 尼 =0. 705, RMSE =5.799 mS + cm， 
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4 3 种 卫星 各 波段 重要 性 


Fig.4 Importance of each band of the three satellites 


表 4 加 入 三 维 光谱 指数 土壤 EC 估算 结果 


Tab.4 Results of soil EC estimation for adding three-dimensional spectral index 


建 模 集 验证 集 
卫星 光谱 形式 参量 
RMSE / mS * em ^! R? RMSE / mS + em ^! R? AIC 
Landsat 8 OLI R 9 4.160 0.757 5.680 0. 750 35.086 
FD 9 3.786 0.813 4.348 0. 835 32.765 
SD 9 3.517 0. 834 5.276 0.777 34.447 
ABS 9 4.075 0.778 5.601 0.725 34.965 
MSC 9 3.645 0. 809 5.459 0. 804 34.743 
Sentinel 2 R 15 3.908 0. 803 5.873 0. 795 47.37] 
FD 15 3.192 0.875 4.817 0. 843 45.656 
SD 15 3.158 0. 865 4.596 0. 857 45.247 
CR 14 4. 186 0. 803 5.799 0. 705 45.267 
ABS 15 3.912 0. 803 5.383 0. 768 46.937 
MSC 15 3.308 0. 888 4.392 0. 834 45.853 
Sentinel 3 R 23 4. 306 0.823 5.986 0. 635 63.543 
FD 23 3.381 0.844 4.636 0. 833 61.324 
SD 23 2.136 0.956 2.986 0.935 57.500 
CR 23 3.271 0.889 5.142 0. 725 62.223 
ABS 23 4. 269 0.801 5. 680 0. 696 63.087 
MSC 23 3.494 0. 859 4.177 0. 840 60.417 
AIC 245.267) ,各 光谱 形式 的 估算 精度 依次 为 SD> 3 讨论 


FD > MSC >R > ABS > CR。Sentinel 3 卫星 基于 SD 
预 处 理 的 土壤 EC 估算 精度 最 佳 (R^ = 0.935, 
RMSE =2.986 mS * cm ,AIC 2 57.500) ,基于 及 的 
土壤 EC 估算 精度 最 低 (R =0. 635, RMSE =5. 986 
mS * em ,AIC =63. 543) ,各 光谱 形式 的 估算 精度 
依次 为 SD > MSC > FD» CR > ABS >R, 综合 比较 


综合 表 3、 表 4 分析 后 认为 ,微分 预 处 理 方式 对 
于 提高 3 种 卫星 土壤 EC 估算 精度 具有 极 好 的 效 
果 ,Sentinel 3 卫星 为 土壤 EC 估算 的 最 佳 选择 ,由 于 
其 27 d 的 重 访 周 期 以 及 300 m 的 空间 分 辨 率 '* , 较 
适用 于 大 尺度 乃至 全 球 范围 内 土壤 盐 渍 化 监测 与 数 


RIR Sentinel 3 卫星 依旧 为 土壤 EC 估算 最 佳 卫星 ， 
其 次 为 Sentinel 2 ,最 后 为 Landsat 8 OLI。 图 $ 为 3 
种 卫星 基于 GBRT 的 土壤 EC 最 佳 估算 结果 图 。 


字 化 制图 ;Sentinel 2 卫星 具有 空间 分 辩 率 高 (最 高 
为 10 m) . 重 访 周期 短 (5 d) 的 优势 1, 对 于 土壤 
盐 渍 化 精细 分 类 以 及 快速 动态 监测 具有 一 定 潜 力 ; 
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图 5 基于 GBRT 的 土壤 EC 验证 效果 


Fig.5 Soil EC verification effect based on GBRT 


Landsat 8 OLI 卫星 虽然 相对 前 两 者 土壤 EC 估算 精 
度 较 低 ,但 最 佳 R 也 达 0. 835 ,加 之 其 数据 连续 性 
较 强 ,可 以 为 区 域 土壤 盐 渍 化 动态 变化 及 预测 提供 
HY HE SCPE 

利用 地 面 实测 VIS-NIR 数据 通过 不 同 卫星 的 波 
谱 响 应 函数 模拟 各 自 宽 波段 数据 ,从 而 构建 土壤 EC 
估算 模型 ,从 理论 上 为 土壤 EC 的 估算 与 土壤 盐 涡 
化 的 监测 提供 了 科学 指导 。 但 是 卫星 实际 观测 过 程 
中 还 会 受到 大 气 以 及 地 物 单元 复杂 等 因素 的 干扰 ， 
造成 卫星 实际 宽 波 段 数 据 与 模拟 数据 存在 差异 ， 
此 模型 的 普 适 性 还 需 大 量 数据 及 实验 进行 验证 , 同 
时 今后 需 增 加 卫星 传感器 类 型 ,并 对 土壤 其 他 属性 
进行 对 比 研究 ,以 期 探寻 最 佳 的 土壤 属性 估算 模型 
及 确定 最 适 卫 星 观测 平台 。 


人 


4 结论 


本 研究 通过 模拟 3 种 卫星 宽 波 段 数 据 , 采 用 5 
种 预 处 理 方式 ,结合 GBRT 算法 构建 艾 比 湖 土壤 EC 
估算 模型 ,并 构建 TDSI 充分 挖掘 波段 间 的 相互 关 
JR ,最 大 程度 上 保留 土壤 EC. 的 敏感 信息 ,有 效 提升 
了 建 模 精度 。 研 究 发 现 ,在 不 同 土壤 EC 条 件 下 ,3 
种 卫星 宽 波段 模拟 数据 具有 相似 的 光谱 趋势 , 均 在 
红 、 近 红外 附近 有 较 高 反射 率 。 构 建 的 TDSI 与 土壤 
EC 相关 性 基本 均 在 0.4 以 上 ,保留 了 与 土壤 EC 高 
度 敏感 的 红 、 绿 、 蓝 、 近 红外 、 短 波 红外 波段 信息 。 
GBRT 算法 对 于 土壤 EC 估算 表现 出 较 高 的 预测 能 
7] ,Landsat 8 OLI Sentinel 2 „Sentinel 3 最 佳 估 算 R 
分 别 为 0. 831,0. 847 ,0. 903 ,在 加 入 TDSI 后 ,精度 
得 到 显著 提升 ,3 种 卫星 最 佳 估算 R 也 分 别提 高 到 
0. 835 .0. 857 .0.935 ,各 卫星 预测 效果 依次 为 Senti- 


nel 3 > Sentinel 2 > Landsat 8 0LI。 本 人 研究 为 土壤 盐 
渍 化 监测 提供 一 种 新 的 思路 及 方法 ,为 针对 不 同 尺 
度 盐 渍 化 灾害 预测 防 控 提 供 科 学 参考 。 
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Abstract; As one of the important environmental problems in the global arid region, soil salinization not only re- 
stricts agricultural production, but also causes serious land degradation. This problem has attracted wide attention 
from many scholars at home and abroad. In the recent stage,the monitoring and evaluation of soil salinization is 
mainly based on the use of different remote sensing data combined with different model inversion methods , and few 
research was seen on the integration of ground measured spectra and satellite remotely sensed data to obtain soil 
salinization information. Therefore ,this study explores the effect of different satellite broadband simulations on soil 
conductivity and reflects soil salinization. In this study ,the Ebinur Lake Wetland Nature Reserve in Xinjiang Uygur 
Autonomous Region, China was taken as an example. 40 soil samples were collected in the field , and the soil con- 
ductivity (EC) and visible-near-infrared ( VIS-NIR) spectral data were measured. The broadband data of the Land- 
sat 8 OLI, Sentinel 2 ,and Sentinel 3 satellites were simulated by spectral response techniques. The three-dimension- 
al spectral index was constructed by simulating wide-band data and its six spectral forms after five pretreatments ; 
first-order differential ( FD ) , second-order differential (SD) , continuum removal (CR ) , absorbance conversion 
( ABS) ,and multivariate scattering correction (MSC). Three satellite soil EC estimation models were established by 
combining the gradient boosting regression tree (GBRT) algorithm, and the accuracy of the model estimation was 
compared after adding the three-dimensional spectral index. Then the soil estimation potential of the three satellites 
was analyzed. The results show that; Landsat 8 OLI, Sentinel 2, and Sentinel 3 have similar spectral trends under 
different soil EC conditions. With the increase of soil EC ,the reflectance of each band of the three satellites is corre- 
spondingly increased ,and both have high reflectance in red and near infrared. The correlation between TDSI and 
soil EC was basically above 0. 4, and Sentinel 3 had the highest correlation based on SD pretreatment, which was 
— 0. 720. TDSI retains the information of red, green, blue , near-infrared and short-wave infrared bands sensitive to 
soil EC , which indicates that TDSI increases the effective information of soil EC while taking into account the remote 
sensing mechanism, which has certain scientific significance. GBRT algorithm has outstanding performance for soil 
EC estimation. Landsat 8 OLI, Sentinel 2, Sentinel 3 based on the estimation model of broadband simulation data 
produced the best prediction accuracy R^ as 0. 831,0. 847 and 0.903 respectively. After adding the TDSI ,the esti- 
mation accuracy of different spectrum forms of the three satellites has been significantly improved. The correspond- 
ing prediction accuracy A^ of the three satellites was increased to 0. 83 5,0. 857 and 0. 935 , respectively. The com- 
prehensive analysis found that Sentinel 3 had the best estimation effect for soil EC ( R^ =0. 935, RMSE = 2. 986 
mS + em ^! , AIC — 57.500) ,followed by Sentinel 2( R? 20. 857 ,RMSE =4.596 mS + cm ^! , AIC =45. 247) ,and fi- 
nally Landsat 8 OLI (R° 20. 835, RMSE =4.348 mS + cm ^! , AIC 2 32. 765). In this study ,the spectral response 
technique combined with GBRT algorithm has been verified that Landsat 8 OLI, Sentinel 2 and Sentinel 3 satellites 
have a good estimation effect on soil EC in arid area. TDSI can deeply mine the synergy information between bands 
and improve the prediction accuracy of the model. The combination of the two methods could provide favorable guid- 
ance for quantitative monitoring and prevention of soil salinization in arid areas. 

Key words: spectral response; soil conductivity ; spectral pretreatment; three-dimensional spectral index; gradi- 


ent boosting regression tree 


